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> 1 INTRODUCTION AND MOTIVATION 



ABSTRACT 

We investigate the potential of the upcoming LOBSTER space observatory (due circa 
2009) to detect soft X-ray flashes from shock breakout in supernovae, primarily from 
Type II events. LOBSTER should discover many SN breakout flashes, although the 
number is sensitive to the uncertain distribution of extragalactic gas columns. X- 
ray data will constrain the radii of their progenitor stars far more tightly than can 
be accomplished with optical observations of the SN light curve. We anticipate the 
appearance of blue supergiant explosions (SN 1987A analogs), which will uncover a 
population of these undcrluminous events. We consider also how the mass, explosion 
energy, and absorbing column can be constrained from X-ray observables alone and 
with the assistance of optically-determined distances. These conclusions are drawn 
using known scaling relations to extrapolate, from previous numerical calculations, 
the LOBSTER response to explosions with a broad range of parameters. We comment 
on a small population of flashes with 0.2 < z < 0.8 that should exist as transient 
background events in XMM, Chandra, and i?05^T' integrations. 

Key words: supernovae: general - X-rays: bursts - stars: fundamental parameters - 
shock waves - instrumentation: detectors. 



A core collapse supernova produces no electromagnetic radi- 
ation until its envelope is completely consumed by the explo- 
sion fireball. This phase ends, however, with a brilliant flash 
of X-ray or extreme ultraviolet photons heralding the arrival 
of the fireball shock at the stellar surface. The "breakout" 
flash is delayed in time, and vastly reduced in energy, rela- 
tive to the neutrino and gravity-wave transients produced by 
core collapse. However, it conveys useful information about 
the explosion and the star that gave rise to it. 

We are motivated to explore the detection and inter- 
pretation of such flashes by several points. First, break- 
out flashes provide strong constraints on some properties 
of their presupernova stars that are poorly determined from 
the optical light curve. Second, whereas previous searches 
for these flashes have failed, the planned LOBSTER experi- 
ment should routinely discover them. Third, this experiment 
will report (within minutes) the location of almost any core- 
collapse supernova (some flashes are extincted) within its 
fleld of view to a distance of many megaparsecs. Data will be 
downlinked to Earth every 90 minutes and analysed via an 
automated process. This allows rapid optical follow-up and, 
in the unlikely case of a very nearby explosion, correlation 



with neutrino and gravity-wave signals. Early warning and 
the precise timing of explosion are valuable for the interpre- 
tation of the optical light curve: for instance, in calibration of 
the Expanding Photospheres Method for distance determi- 
nation. Fourth, flash-selected surveys of Type II supernovae 
are less biased by competition with light from galactic nuclei 
than are optical surveys. Finally, there exists the possibility 
that LOBSTER or a subsequent instrument could detect the 
signature of asymmetry in the supernova explosion through 
its effects on the breakout flash. 
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In this paper we consider primarily what can be de- 
termined about the stellar progenitor from the X-ray flash 
alone, from the flash and an optical distance, or from the 
delay separating the flash from the core collapse. Our in- 
vestigation is based on previous numeri cal calculations of 
the br e akout flash spectrum, prima r ily bvlKlein fc ChevalieJ 
fl978VEnsman & Burrowi^ (1199211 . iBlinnikov et al.l Jl99d) ' 
and Blinnikov ot al. (200QJ, and on the analytical approx i- 
mations and scalings derived bv lMatzner fc McKed lll999t) . 
Our consideration of observational points relies on technical 
details kindly provided by Nigel Bannister and the LOB- 
STER science team. 
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1.1 Limitations of Optical Light Curves for 
Constraining SN Properties 

Supernova properties provide a datum at the end of stel- 
lar evolution. Collectively, such data from many stars re- 
veal how the intrinsic properties (e.g., initial mass, rotation, 
metallicity) and extrinsic properties (e.g., binary member- 
ship and separation) of massive stars affect their evolution. 
It is currently uncertain, for instance, how many super- 
novae come from relatively compact blue supergiant pro- 



genitors and therefore produce underluminous l i ght curves 
' illl988t iFiUipenkc 



like that of SN 1987A llSchmitz fc Gaske: 



1988|: _ lYoung fc Br anch' 'l98Sl: Ivan den Bergh fc McClure 



19891: IScha,efetl ll 996 V This arises both from a bias against 
detecting small progenitors, and a bias to overestimate the 
radius, as discussed below. X-ray observations have the po- 
tential to mitigate the former bias (Table |1J and eliminate 
the latter 

U s ing the scaling relations of iLitvinova fc NadezhirJ 
|l98l,|Hamu3 i2003l) estimated progenitor properties from 
the plateau properties of the 13 best-measured Type IIP SNe 
(the only ones for which sufficient observations exist; prior to 
his work, only SN 1969C had been analyzed in this fashion). 
Based on measurement error, he assigned la uncertainties 
of roughly 0.26 dex (a factor of 1.8) to the derived radii (and 
similar uncertainties to mass and explosion energy). Several 
additional sources of uncertainty degrade the determination 
of progenitor radius derived from optical data. 

For one, there is an uncertainty in the underlying dis- 
tr ibution of ejecta density and o pacity. The scaling relations 
of ILitvinova fc NadezliinI lll985ll are derived from a suite of 
numerical simulations involving an explosion in a simpli- 
fied version of a red supergiant (RSG) model, in which the 
progenitor density scales with radius as po oc (1 — r/R-i)^^^ , 
where is the size of the star, f ODOv (1993), m comparison, 
developed analytical formulae for plateau-type light curves, 
assumi ng a constant density d istribution in the ejecta. A\- 
though [Litvinova fc NadezhirJ ' s results are nominally more 
precise, being derived from a hydrodynamical calculation of 
a model star, neither set of formulae account in detail for 
the presupernova structure of the progenitor - which de- 
pends, for instance, on the mass ratio between the star's 
mantle and its hydrogen envelope. (Simple formulae for the 
ejecta distributions from supergiants have been presented 
by iMat zner fc McK ee (1999), but have not yet been used 
to predict light curves). For a typical plateau light curve 
(absolute V magnitude peaking at -17 .5""^^ at 70 days; 
early photospheric ve locit y 7000 km/s; IPopovI 119931) . the 
ILitvinova fc NadezhirJ and lPopovl determinations of 7?* dif- 
fer by 0.28 dex or a factor of 1.9, which is comparable to 
the observational uncertainties. This serves as an estimate 
(although probably an overestimate) for the errors due to 
an unknown ejecta distribution. 

Additionally, there is a bias toward large values of _R* 
due to the contamination of the plateau peak by the ra- 
dioactive d«ca^;_£f^2Co_2rodAicedm_the e xplosion. The for- 
mulae o f iLitvinova fc NadezhinI and lPopovl ^9^) as- 



sume the light curve to be powered by heat deposited in 
the ejecta by the explosion shock itself. At the time that 
photons can diffuse out of the ejecta, this heat supply has 
been adiabatically degraded during expansion; the fraction 
remaining is inversely proportional to the initial stellar ra- 



dius. In contrast, Co, with a half-life of 77.3 days, releases 
its energy during the plateau. If the initial radius is small 
enough, cobalt decay will compete with or even dominate 
shock-deposited heat. To detect this effect, one must monitor 
the late-time decay of the light curve to estimate the cobalt 
contribution. Nevertheless, the derived progenitor radii are 
easily corrupted, because there is no simple way to sub- 
tract the radioactive contribution from the light curve; nor 
is this possible, if the ^^ Co lum i nosity exceeds that from 
shock heating. Using the IPopovI Jl993h formulae, we find 
that radioactivity dominates when 
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where R = 5OOR500RQ (an alternate, R = 5OR50RQ, will 
be used later for blue supergiants); the ejecta mass is Mej = 
lOAfoj^ioAfQ; the ejecta opacity is k = 0.34ko.34 cm'^ g~^; 
and the explosion energy is Ein = 10^^ £51 erg. 

For the fiducial red supergiant (i?5oo = -E51 = Afio = 
K0.34 = 1), equation Q indicates that O.21Af0 or more of 
^^Co would be required to dominate the plateau l uminos- 
ity. (This is a large but not impossible amount: see iHamuvl 
[2003] 's Table 4). In contrast, reducing the initial radius to 
50i?Q, only 0.038 Mq of ^^Co is needed to dominate the 
luminosity. Clearly, one cannot detect the existence of blue 
supergiant (BSC) progenitors from light curve plateaus if 
there is even a minute amount of radioactive cobalt. This 
can also be seen by solving equation Q for the value of R 
below which ^"Co dominates: 
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This local, power-law solution to the implicit equation Q 
is reasonably accurate for typical doses of ^®Co. The scaling 
Afseco c< RI "^^ indicates that ^®Co contamination is an issue 
even for red supergiants. 

These difficulties prohibit the use of a catalogue of 
plateau-phase SN light curves to derive the population of 
red versus blue supergiant progenitors. Further, such a cat- 
alogue would necessarily be biased by the fact that larger 
stars and those with more ^^Co have intrinsically brighter 
optical displays. 

There are other means for placing constraints on J?* 
optically. One possibility is to locate the progenitor star 
in previous optical observations and to estimate i?* from 
the observed luminosity and spectral type (see, for exam- 
ple, ;van Dvk ct al. 2003). Of course, this method hinges 
on whether the progenitor identity can be established from 
archival data. 

Alternatively, sufficiently early observations may catch 
the SN light curve in the phase of self-similar spherical dif- 
fusion, during which the light curve is powered by shock- 
deposited heat escaping from the highest- velocity ejecta 
(|Che^iie3[l993). This theory constrains the combination 
Ri, E^^^^ M~j^'*'' rather than i?* alone. For instance, an ap- 
plication to_tlie_ea£l^_observations of the Type Ib/c SN 
1999ex JStritzinger et alj I2OO2I) yields an estimate ii* ~ 
7.0E^i^^M~yj^°RQ. To yield R^,, this method requires a 
knowledge of the time of explosion, observations of the SN 
in its first days, and estimates of the distance, explosion 
energy, and ejected mass. 
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As we shall show, is the parameter most tightly 
constrained by the breakout flash. X-ray observations are 
thus independent of and complementary to opti cal follow- 
up, a nd provide the early warning needed for IChevalierl 
il992h 's method. Indeed, X-ray and optical constraints can 
be merged to provide a more complete picture of the discov- 
ered SNe. 



1.2 Previous work 



Shock breakout flashes were predicted bv lColeat^ il96d) as 
a so urce for (the then undete cted) 7-ray bursts. 

iKlein fc ChevaheJ il978l) carried out radiation hydro- 
dynamical calculations for the explosions of red giant stars, 
predicting breakout flashes detectable by the soft X-ray tele- 
scope HEAO-1. Unfortunately, the field of view of HEAO-1 
was insufficient to detect any breakout fiashes in the limited 
duration of the experiment (Klein ct al]|l979ll . 

The explosion of SN 1987A, and the realization that 
archival images of its progenitor indicated a blue rather than 
a red supergiant sta r, stimulated a reana l ysis o f supernova 



breakout flashes by Ensman & Burrowi 



1 099?) 



and, most 



recent ly, by [Blinnikov et al. ( 1998) and Blinni kov et alJ 
i2000l) . These studies represent an increase in sophistica- 
tion toward the full numerical treatment of this complicated, 
radiation-hydrodynamic problem. They are consistent with 
the observed ionizati on state of circumstcUar gas around the 
SN 1987A remnant C Ensman fc Burrows 1992). 

The physics of shock breakout and the accompa- 
nied flash have also been treated more generally with 
analytical approximations for many years. The hydro- 
dynamical problem - that of a shock wave accelerat- 
ing down a de nsity ramp - was solved in a self-similar 



idealizat i on by iGandel' man fc Frank-K amenetskvl il95d) . 
Sakurail (ll96ol) a n d iGrover fc Hardvl lll96d) aiid used by 
■Colgate fc White! Jl96d) m their study of supernova hy- 
drodynamics. AfteP colgatd il968l) 's initial suggestion, the 
ph ysics of the breakout fla s h was reex amined analytically 
bv llinshennik fc Nadezhin I (Il989l) and iMatzner fc McKed 
Jl999h : see O 

We wish to tie together the numerical and analytical 
work on breakout flashes to predict the detectability by 
LOBSTER of SNe of a wide variety of characteristics. This 
requires a method of scaling the results of a numerical simu- 
lation to a different progenitor mass, radius, or explosion en- 
ergy (or, less important ly, envelope structur e). This is made 
possible by the work of lMatzner fc McKed (1999; hereafter 
MM 99), who developed a general analytical approximation 
that accurately predicts the speed of the explosion shock 
front. Evaluated at the stellar surface, and combined with 
llnislicimik fc Nadezhin (1989)'s theory, this gives an esti- 
mate of the onset, duration, luminosity, and spectrum in 
the flash. Although numerical studies predict the spectrum, 
for instance, more accurately, MM 99 's formulae provide the 
scaling laws required to generalize them. 



1.3 This Work 

In this paper we examine the ability of sensitive, wide- 
field X-ray detectors (in particular, the forthcoming LOB- 
STER instrument) to constrain the properties of supernovae 



through observation of the shock breakout flash. We begin, 
m 31 with a brief overview of the physical processes which 
generate the breakout flash. Section |H] describes, in detail, 
how the intrinsic properties of a supernova and its progeni- 
tor star can be gleaned from the characteristics of the X-ray 
flash. We attempt to account for the distribution of Galactic 
and extragalactic interstellar X-ray absorption by employ- 
ing Cappellaro et al. (1997)'s estimate of the distribution of 
optical extinctions toward supernovae; however this distri- 
bution is quite uncertain. A brief discussion of the shock 
travel time and its usefulness in constraining parameters 
is included (2|l; unfortunately, LOBSTER is very unlikely 
to observe a supernova near enough for gravity waves and 
neutrinos to be detected. We include a short discussion on 
the use of instruments other than LOBSTER for detecting 
breakout flashes {^^; we conclude that there may be un- 
noticed breakout flash detections in the archived images of 
the XMM, Chandra, and i? 05^ T X-ray telescopes. In ^ we 
summarize how shock breakout observation can be used as 
a powerful tool for understanding the end stages of stellar 
evolution and underline some caveats. In the Appendix, we 
discuss, in greater detail, the signiflcance of the outer den- 
sity distribution of the progenitor and show how it can be 
characterized by other stellar parameters. We model stel- 
lar luminosity as a function of stellar mass for presupernova 
stars and use this relation to estimate the outer density co- 
efficients for blue supergiants. 

Our analysis improves upon previous work by dealing 
more generally with breakout fiash properties and by mod- 
elling more specifically the observation of bursts. Previous 
work has concentrated on the prediction of the breakout 
fiash from a single star; the method of scaling used in this 
work allows for us to describe fiashes from a broad range 
of progenitors (albeit in less detail). We concentrate on the 
observation of these events by LOBSTER and the recon- 
struction of the supernova properties from the data. 



2 PHYSICS OF BREAKOUT FLASHES 

We briefiy review here some of the physics relevant to 
breakout flashes here; howev er, we refer the reader to 
llmshennik fc Nadezhin | il989l) and MM 99, and references 
therein for a more thorough discussion. 

MM 99's theory for the speed of the explosion shock 
front (their Eq. [19]) reads 



Vs{r) = 0.794 



" Sin ■ 


1/2 


m(r) 


m(r) 




r-3po(r) 



(3) 



where _Bin is the explosion energy, r labels position within 
the progenitor, and p(i(r) and m(r) are the progenitor's den- 
sity and ejected mass interior to r. This formula accounts in 
a simple manner for two effects: deceleration as the shock 
sweeps up new material (the flrst bracket), and acceler- 
ation in regions of sharply declining density (the second 
bracket). The prefactor is matched to a self-similar solution; 
^an, Matzncr, fc McKee ( 2001i) have given a slightly more 
accurate form in which this coefficient depends on the shape 
of the density distribution. Likewise, the exponent 0.19 on 
the second term can be adjusted (very slightly) depending 
on the slope of the density profile near the stellar surface. 
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Nevertheless, equation is accurate within 10%, often 3%, 
as it stands. 

The breakout flash itself occurs because supernova 
shocks are regions of radiation diffusion, and when they 
get too close to the stellar surface, photons diffuse out into 
space. Such shocks are driven predominantly by radiation 
rather than gas pressure, and this is increasingly true as 
they accelerate through the diffuse subsurface layers. The 
shock thickness (Is) is therefore determined by the require- 
ment that photons must diffuse upstream across the shock 
(in a time ~ Tsh/c, where Ts is the optical depth across Is) 
as fast as fluid moves downstream across it (in a time h/vs). 
Equating these. 



fVs 



(4) 



We have included a factor / ~ 1 — 3 which is uncertain at 
this level of approximation because the shock is a smooth 
structure; MM 99 assume f ~ 1, arguing its effect is unim- 
portant. We retain / to show its effect, but note that it 
always appears below as a coefficient modifying the X-ray 
opacity k. 

The quantity Ts decreases as the shock accelerates, but 
the optical depth to the surface (r) decreases more rapidly. 
Once T < Ts, there is insufficient material to contain the 
shock, and the photons that constitute it will escape and 
produce the flash. 

To specify the flash's properties, one requires a formula 
for the subsurface density profile; MM 99 motivate 



(5) 



where n is related to the polytropic index 7p = dlnp/dlnp 
via "fp — 1 + 1/n. Hence n ~ 3/2 in an efficiently convecting 
atmosphere (i.e., in red supergiants) and n ~ 3 in a radiative 
atmosphere of uniform opacity (blue supergiants). Evaluat- 
ing equation JSJ at r = R*/2 results in pi = p(_R*/2); thus, 
pi represents the half-radius density. The outer density dis- 
tribution of a star can be described by the ratio pi/p*, where 
Pi, is the characteristic density Mcj/R^ (see i^for more de- 
tail). 

Finally, one requires the relevant opacity k in the ex- 
pression T = J pK, dr. Since the post-shock temperature is 
~ lO'^ - 10® K, MM 99 assumed k would be dominated 
by electron scattering — hence, k ~ 0.34 cm^ g~^ for so- 
lar metallicity. We have investigated this assumption for the 
relevant densities and te mperatures using the OPAL tables 
ilglesias fc Roger3ll99(J) . and find it perfectly valid. 

The results of setting Ts — t are given by MM 99 in their 
§5.3.3; we will not repeat them here. Their formulae (36), 
(37), and (38) give the characteristic temperature Tsc (hence, 
photon energy ~ 2.7kBTsc), flash energy Esc, and radiation 
diffusion time tsc, respectively (the latter sets a lower limit 
for the observed flash duration). We do, however, wish to 
note several important features of these formulae (Table 0. 

First, these quantities depend most strongly on the pro- 
genitor radius _R*, then, decreasingly, on Ein, the total mass 
of ejecta (Mej), k, and flnally, quite weakly on the outer den- 
sity coefficient pi/p*. For this reason, we will ignore varia- 
tions of the structural parameter pi /p* in the current paper, 
and instead concentrate our attention on Ri,, the parameter 
to which flashes are most sensitive. We caution, however, 
that the effect of pi /p^, should be considered in cases where 



the outer density profile can be determined. An example 
would be a radiative envelope of known lu minosity; as dis- 
cussed bv iTan. Matzner. fc McKed i200 j) . pi/p* becomes 
very small if the star's luminosity approaches its Eddington 
limit. In the Appendix we present alternative breakout scal- 
ing relations that incorporate analytical approximations for 

Second, there is a correlation between the photon tem- 
perature Tse and diffusion time tsc induced by the physics 
of shock breakout. Since photon pressure dominates, Tsc 
must satisfy aTst/3 = (6/7)pff at the point of breakout. 
And, since the diffusion time of the shock equals its cross- 
ing time, tsc = Is/ Vs. These expressions, combined with 
r = J Kp dr — Ispn/iji + 1) and r = c/{fvs), give 



18(n + l)c 
7 fan 



(6) 



As we have noted, k ~ 0.34 cm^ g~^ for envelopes of roughly 
solar composition. Insofar as tsc and Tse can be measured, 
equation Q provides a means to determine observationally 
whether the duration of an observed flash is dominated by 
diffusion across the shock front. 

Alternatively, the flash duration may be dominated by 
the finite difference in light tra vel time between the centr e 
and the limb of the stellar disk jEnsman fc Burrow jll99'3l . 
We estimate this as 

(7) 

although one should keep in mind that the emission pattern 
at each point on the stellar surface (i.e., limb darkening) 
will make the effective value of tc slightly smaller. In gen- 
eral, the observed flash is a convolution of the diffusion and 
propagation profiles; hence the total duration is 



^ — (^sc ^c) ; 



(8) 



although neither profile is actually Gaussian. 

Returning to equation Q, one sees that k will be un- 
derestimated if one substitutes t for tsc, provided tc > t. 
Likewise, the characteristic photon temperature can be in- 
creased relative to Tse by selective absorption along the line 
of sight. If one derives k <C 0.34 by substituting observed 
quantities in Eq. Q, then either tc > tso, or extinction is 
significant, or (quite likely) both. Unfortunately, LOBSTER 
will not view bursts long enough ( ^13.111 to apply this test in 
the cases where it would be useful. We concentrate below on 
other means for constraining SN progenitors. 



3 SUPERNOVA PARAMETERS FROM X-RAY 
CONSTRAINTS 

There are five explosion parameters one would wish to know: 
the intrinsic parameters J?*, _Ein, and Moj; the distance D; 
and the obscuring column A^h. In addition, the structural 
parameter pi/p* plays a minor role. As we shall see in the 
Appendix, pi/p* is determined approximately by the stel- 
lar luminosity (in BSGs) or the mass fraction in the outer 
envelope (in RSGs). In each case, however, pi/p* is roughly 
constant at a characteristic value. 

As described in 321 the intrinsic parameters of the explo- 
sion determine the characteristic properties of the breakout 
flash: its colour temperature Tsc, its total energy Esc, and its 



SN Properties from Shock Breakout 5 



Table 1. Scaling equations from MM 99. These should be read as exponents in formulae, e.g., tge = lO"'^* [/k/(0.34 cm^ g ^)] 



cgs units. 
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-0.086 
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0.21 
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Figure 1. Construction of the LOBSTER instrument response 
to a blue supergiant SN at a distance of 2Dniin = 50 Mpc. Eight 
events are expected at this distanc e (Eg. P]*). Top p anel: in- 
terstel lar cross section per H atom iWilms et al.l l2000Y Second 
panel: iBlinnikov et al.l j^OO)'s model for the breakout flash of 
SN 1987A, scaled with MM 99's scaling relations (solid) and ex- 
tincted by a typical column (dashed). Third panel: LOBSTER ef- 
fective area. Bottom panel: counts for the unextincted (solid) and 
extincted (dashed) flash, binned in 20% bins (the energy uncer- 
tainty of the detector). Also shown are the three energy channels 
we use to characterize colour. 



duration t. In theory, the full breakout spectrum provides a 
large number of parameters. However, we take the view that 
only these three are independent. If so, then it is not pos- 
sible to reconstruct J?*, Ein, and unless pi /pi, is held 
constant. Variations of pi/p* around its characteristic value 
leads to some uncertainty in the other parameters. 

The X-ray observables depend on D and A^h in addi- 
tion to t, Tsa, and Eso- Since there is negligible redshift at 
the distances in question, D affects only the X-ray fluence 
(oc Ese/D^ for given Tae). For this reason, E^c cannot be 
known unless D can be ascertained from optical foUowup 
observations. 

Given the ~ 20% uncer tainty in L OBSTER'S deter- 
mination of photon energies jBannisteJliooA . it provides 
roughly nineteen independent energy channels ( tl3.11 : with 
the burst duration, this makes twenty observables per flash. 
However, we use only four parameters in our analysis: three 
channels (defined in i|3.2L plus the burst duration t. 

The motivation for this choice is twofold. First, there 
are five desired intrinsic parameters, two of which {Esc and 
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Figure 2. As in Figure ITl but fo r a red s upergiant model. The 
breakout flash was calculated by IBlinnikov et al.. a99&) for SN 
1993J; see 



D) cannot be disentangled without optical foUowup; there- 
fore only four X-ray observables need be used. Second, if 
supernovae occur homogeneously within the detection vol- 
ume, the number of bursts observed scales as the limiting 
fiuence according to A'' oc(fiuence)~^'^^. The median burst 
fiuence is therefore only 60% above the detection threshold, 
and cannot be described in fine detail. 

After describing the LOBSTER instrument (EHJ. we 
describe our method for simulating X-ray observables ( tl3.21 
and our strategy for constraining supernova parameters with 
these four observables (S H3.3t03l . 



3.1 Model of the LOBSTER instrument response 

LOBSTER is a proposed wide-field soft X-ray detector, cur- 
rently envisioned to be deployed on the International Space 
Station in 2009. LOBSTER has a field of view of 162x22.6 
degrees and would continually scan the sky as it orbits the 
Earth every 90 minutes. Doing this, it would image ^ 95% 
of the sky once every orbit. LOBSTER has a spatial res- 
olution of approximately 4 arcmin, and is sensitive to soft 
X-rays of energy 0.1 - 3.5 koV from sources with fluxes of or- 
der 10~^^ erg cm~'^ s~^ (.Priedhor skv ct al. 1996). Incident 
photon energies are recorded with an accuracy of roughly 
20%, implying ~ 19 effective spectral channels. 

LOBSTER focuses X-rays through grazing angle reflec- 
tions by an array of microchannels. This arrangement pro- 
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Table 2. SN detection distance and rate versus model type. The factor f^hs represents the unknown popu- 
lation of each subclass (Eq. In each model, Moj = 15Mq, = 10^^ erg, and re = 0.34 cm^ g~^, unless 
otherwise noted. 



Progenitor 






^max 


Detections 


Note 


Type 


(1021 cm-2) 


iR&) 


(Mpc) 


per year 




BSG 


1 


50 


260 


- 320/obB 


SN 1987A analogue 


BSG 


5 


50 


86 




SN 1987A analogue 


RSG 


1 


500 


420 


~ 1400/„bB 




RSG 


5 


500 


66 






RSG 


1 


500 


800 


~ 10000/obs 


SN 1993J analogue; M^j = 2Mq 


Wolf-Rayet core 


1 


5 


58 


~ 1/obs 


BSG scaled to SN Type Ib/c analogue; 












Mcj = 3Mq, k = 0.2 cm^ g-^, pi/p* = 0.1 



duces a cruciform focal spot on the detector, with a bright 
central focus spot at the intersection of two dimmer, orthog- 
onal arms; there is also a diffuse background of unfocused 
photons. The instrument's effective area is a complicated 
function of incident photon energy, being limited at low en- 
ergies by detector window absorption and at higher energies 
by decreased reflectivity. Shown in the third panel of Figure 
Qis the effective area for the central focus, from data pro- 
vided by Nigel Banni ster (2003, private commu nication). For 
a full discussion, see IPriedhorskv et al.Hl99d) . 

The instrument's field of view should allow for success- 
ful detection of numerous shock breakouts. One can define 
the instantaneous volume of view 

Vvicw = ^^-^max 

where -Dmax is the maximum distance at which a certain 
flash can be detected. The field of view, Q, is a narrow swath 
of the celestial sphere: Q = 2(162°) sin(22.6°/2) = 1.11 sr. 

The instantaneous volume of view is appropriate for 
events that are shorter than the dwell time fdwcii of the ex- 
periment. The dwell time for LOBSTER has a minimum 
value of roughly (90 min)(22.6° /360°) = 340 s, with longer 
times occurring near the poles of the orbit, to a maximum of 
tdweii = 1800 s jSannistcr 2003). We will adopt the typical 
value of 400 s for tdwcii- Longer events can be found even 
if their beginnings are not observed, leading to a larger ef- 
fective solid angle Qctt — (1 + t/tdwcu)^, until the event is 
viewed in more than one scan of the sky (although this is 
not expected for breakout flashes) . This increase in sensitiv- 
ity for long events comes at a cost: the burst duration t is 
no longer directl y observable, and c a n only be constrained. 

Note that iKlein fc ChevalieJ lll978l) assumed that 
breakout flashes would outlast the 30 s dwell time of HE AO 
1, and calculated their detection rate accordingly. LOB- 
STER'S dwell time, in contrast, is an order of magnitude 
longer; only RSGs (with > ITORq) have t > tdwcii for 
this experiment. 

ICappellaro et all 1^9^ quote a rate of 0.71h^ Type 
II supernovae per century per 10^" solar luminosities in 
the B band. There a re 2 x 10^hLp,Mpc~ ^ in th e B band 
JPukugita et alJll998h : from ISoergel et al.l J2003l) we adopt 
h = 0.71 (the Hubble constant in units lOOkms"^ Mpc~^). 

If a fraction fobs of Type II supernovae belong to a sub- 
class of interest, and if this subclass is observable to distance 
Dnia.x, then the number of flashes of this type expected in 
the total integration time of the instrument's life, tohs, is 



/o 



D 

mm {tabs 



(9) 



The quantity -Dmin(tobs), a function of the duration of the 
experiment, is the minimum distance at which one would 
expect, statistically, to observe one SN II of any progenitor 
type: 



-Dmiii(tobs) ~ 26 



3 yr / tdwcii 

^obs V i ^dwcll 



1/3 



Mpc. 



(10) 



Although this is comparable to the distance to the Virgo 
cluster, the universe is reasonably homogeneous on scales 
D > -Dmin. (The local overdensity may reduce the actual 
minimum distance by ~ 25% relative to Dmin; however we 
neglect this effect in its definition.) 

For the LOBSTER mission, tobs ~ 3 years; Table H 
shows -Dmax aud A^obs for various progenitor types, including 
the canonical SNe 11 (a RSG progenitor with ALj — 15Mq , 
Ri, = 5OOJ?0, Ein = 10^^ erg, and a BSG progenitor with 
Moj = 15Mq, R^, = 50i?Q, Ein = 10^^ erg). Although the 
number of events observed depends strongly on parameters, 
LOBSTER should observe between ~ 10^ and 10^ SNe II 
and of order one SN Ib/c over a three year lifespan. In tl3.2.1l 
we refine these estimates by averaging over the distribution 
of absorbing columns A^h; see Figure |3 and Table |1| 

Note that in Table |5| we have included Type Ib/c su- 
per novae, which occur five times less frequently than Type 
lis iCappellaro et aP Il999l) and hence have Dmax = 44 
Mpc. Their observability is rather strongly dependent on 
the rather poorly known initial radius, ejected mass, and 
explosion energy: for these events, A'obs oc (Ese/Tsc) 



3/2 



^iii ' R-i, (quite approximately, because Aaf{{hv) is not 
flat and extinction is not completely negligible). In the table, 
we denote by /obs an unknown fraction with the properties 
given in Table |5| 



3.2 Prediction of Breakout Flash Properties 

We describe in this section the prediction of X-ray observ- 
ables from the properties of a model explosion. We begin 
by calculating the intrinsic parameters Esc, Tsc, and tsc us- 
ing the equations presented by MM 99 as shown in Table 
We choose from the literature a numerical calculation of the 
breakout spectrum as close as possible to the desired explo- 
sion, then shift the spectrum to make its total energy and 
mean photon energy conform to these model predictions. 
Unfortunately, there are very few published calculations 
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Table 3. Parameters used in scaling breakout spectra. 



Progenitor 


Mo, 


(Mq) 


R4Rq) 


Pi/p* 


type 










BSG 


1 


- 35 


20 - 100 


0.2 


RSG 


1 


- 35 


50 - 1000 


0.5 



of brea kout flashes. For red supergiants there is the calcula- 
tion bv lKlein fc Che^^^r 1|197&') , and t he more recent mode l 
for SN 1993 J jBlinnikov et"anil998t) . iKlein fc ChevalieJ 's 
model incl udes a hard tail of higher-energy photons that 
is absent in lBlinnikov et eiLI 's model, presumably due to the 
more so phisticated radia t ion tr ansfer employed in the latter. 
Because iBlinnikov et alJ (^9^ employ multi-group transfer 
calculations, we adopt this calculation as the fiducial RSG 
flash. However, there is insufficient information about the 
progenitor in the literature to calculate Tsc for it using the 
MM 99 equations. When using it, therefore, we enforce that 
the mean photon energy agrees with the value 2.7fcTsc ap- 
propriate for a blackbody of that colour temperature. 

For blue supergiant s t here are two pot e ntial s ources 
of breakout spectra : lEn sman & Burrows! 1^93) and 
IBlinnikov et al.l ll200d) . both of which were calculated for 
SN 1987A. We use the latter, again because the multigroup 
radiation transfer method employed therein is the more so- 
phisticated. In this case the MM 99 formulae could be ap- 
plied to the progenitor model, and predict, within 20%, the 
mean photon energy. We nevertheless enforce a mean en- 
ergy of 2.7kTse when constructing model BSG flashes, to be 
consistent with the RSG case. 

For all Type II supernovae we fixed /k at 0.34 cm^ g~^ 

(©• 

Our three spectral channels were defined as follows: 0.10 

- 0.33 keV (channel 1), 0.33 - 0.54 keV (channel 2), and 0.54 

- 3.5 keV (channel 3). The chosen energy range matches the 
response of the LOBSTER instrument, 0.1 - 3.5 keV ( JO 
Figs.Qand|2l . The lower bound of channel 2 was chosen such 
that all hashes had sufficient channel 1 photons at D = Dmin 
after suffering the expected amount of interstellar extinction 
(A''h ~ 10^^ cm~^); the upper bound was set so that all 
RSG flashes had at least some channel 3 photons prior to 
extinction. 

If Ci,2,3 are the counts in the three channels, we may 
define colour parameters C2-1 = logjQ(C2/Ci) and C3-2 = 
logj^Q(C3/C2). In the absence of an independent distance de- 
termination (see the introduction to these two colour 
parameters and the flash duration t are the only constraints 
on the explosion itself. Unlike t, C2-1 and are affected 

by absorption. 

Figures Q and 121 fllustrate the prediction of the LOB- 
STER instrument response to typical blue and red super- 
giant explosions, respectively. 

3.2.1 Interstellar X-Ray Absorption 

In the X-ray, bound-free absorption by the interstellar 
medium (ISM) is greatest for the softest photons. This leads 
to a "bluening" of the X-ray spectrum as opposed to red- 
dening in the optical. To account for t he efi^ects of absorp - 
tion, we apply the X-ray opacities of I Wilms et al.l (I2OOOI) 




Figure 3. Maximum distance of observability, -Dmax, versus ab- 
sorbing column Afjj for several progenitor models. Models with 
harder spectra are less sensitive to extinction. Also plotted are 
the distances -Dmin within which one supernova is expected in 
three years; one expects ~ D'L^^ / T>^;„ events. The distribution 
of co lumns due to the Milky Way ISchleg el^ Finkbeiner. &: Davisl 
Il998^ and inferred for host spiral galaxies lCaPDeUar^^^ a^^^oW 
are shown (dotted lines: 10th to 90th percentiles; solid dashes: 
25th, 50th, and 75th percentiles). Note that the host column dis- 
tribution is very uncertain. 

and consider a range of hydrogen column densities A^h . 
Milky Way columns JSchleeel. Finkbeiner. fc David [l99l) 
can roughly be described with a log-normal distribution: 

logio(AH/cm-2) ~ 1020-72±0.62 ^^^^y 

Host galaxy column densities toward Type II and Type 
Ib/c supernovae are quite uncertain; however there is evi- 
dence they are often significantl y large r than typical Milky 
Way columns. ICappellaro et all (1199711 estimate blue mag- 
nitudes of extinction toward SNe of both types by searching 
for trends with galaxy inclination (after accounting for the 
type of the host and adopting a fiducial value for the face- 
on extinction). Adopting the model plotted in their Figure 
4 leads^ to a range of Ah from 1.5 x 10^^ cm"^ to 1.3 x 10^^ 
cm~^, in a distribution strongly skewed toward higher val- 
ues. 

These are comparable to the columns through Galactic 
molecular clouds, which may reflect the interstellar context 
in which they explode. On the other hand, molecular gas is 
seen to disperse from OB associations on timescales shorter 
than the typical presupern ova lifetime (lBlaauwll99 J) . More- 
over, ^^^^^Berg^ 1I1991I) argues that chimneys around OB 
associations produce an observed bias toward face-on galax- 
ies. In these two ways, energetic feedback from OB stars may 
reduce Ah in a fraction of the supernova population. 

The Galactic and extragalactic column distributions are 
shown in Figurej^ Since Aobs — -Dmax/^mim if is quite clear 

^ We quote columns in H atoms cm~^, but these are derived from 
magnitudes of extinction in B, and applied to X-ray absorption. 
The end result should therefore be roughly independent of host 
galaxy metallicity. 
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Table 4. SN detection rate averaged over a distribution of extinctions. SN models from Tablel^ the Galactic 
column density distribution I Schlcgel. Finkbcincr. & Davis 1998) and an estimate of the extragalactic column 
distribution (Fig. 4 of ICappellaro et al.lll997ll are used to find the typical distance ("Eg. lilt observable 
number, and column density (averaged over observable flashes, Eg. 1121 toward each type; see Fig.|21 These 
are sensitive to the small-A'^jj tail of the extragalactic column distribution, which we reiterate is quite 
uncertain; see Figure 1^ 



Type 






Detections 


(logioArH(cm-2)^ 




(%) 


(Mpc) 


per year 




BSG (87A analogue) 


50 


96 


17/obs 


21.40 


RSG 


500 


113 


27/obs 


21.32 


RSG (93J analogue) 


500 


223 


210/obs 


21.33 


W-R core (Ib/c) 


5 


37 


0.2/obs 


21.68 



that the statistics of the LOBSTER breakout-flash catalogue 
will be sensitive to the Iow-A'^h end of the host-galaxy col- 
umn distribution. This is demonstrated in Table 0] which 
gives the effective value of Dmax averaged over the extinc- 
tion distribution: 



(-^max) — 



(11) 



where V{< Nn) is the probability that the column den- 
sity is less than Nn- The expected number is then A'^obs = 

/obs(Dmax)/Anin aS iu Eq. © . 

Table El also gives the log- normal value of A''h, averaged 
over the observed bursts: 



(logio Nh) 



(12) 



Although these numbers are subject to our uncertainty of 
V{< Nu), they give a sense of what may realistically be 
observed in a survey. Both estimates strongly weight the 
low end of the A'h distribution, causing the 'typical' column 
to be only ~ 2 - 2.5 x 10^^ cm"^. 

The stellar wind of the progenitor star can have 
a column far in excess of the interstellar values: for 
a mass loss rate ~ 10~^MQ/yr and a terminal ve- 
locity of order the escape velocity, it is ~ 2 x 
W'^^WMQ/M^y/^{50RQ/R^,y/^cm-^. However, this wind 
is fu lly ionized by a sm all fraction of the breakout photons 
(see iLundavist fc Fransson .1988) and cannot significantly 
affect the X-ray spectrum. 

The closest supernova expected in the LOBSTER cat- 
alogue is located Dmin ~ 26 Mpc away, as defined in 
equation assuming three years of observation (RSGs 
with t > tdwoii should be detected at a shorter distance; 
>13.H . If the total number of counts detected from a su- 
pernova at -Dmin is less than 10 (a fiducial number), we 
consider the supernova unobservable. (Unobservable fiashes 
have -Dmax < Dmin-) Likewise, we consider the colour param- 
eters C2-1 and Cs~2 uncharacterized if there are not enough 
counts (< 5 in a necessary channel) to construct them given 
a SN at Dmin- These definitions are employed in Figures 0| 

and|7|to identify excessive values of A^h- 

When the effective area of the LOBSTER instrument is 
taken into account, analysis reveals that all BSG and RSG 
progenitor models with Ein ~ 10^^ erg are unobservable 
when A'^H > 10^^ cm~^. The colour parameter C2-1 is un- 



characterized for all models when A^h > 2 x 10^^ cm ^ and 
6*3-2 is uncharacterized when A^h > 8 x 10^^ cm~^. 



3.3 Constraints From Timing Alone 

It is possible to place constraints on the progenitor radius 
simply by measuring the duration of the shock breakout 
burst and comparing the light travel time to the diffusion 
time for the given range of parameters in Table |3 

One can define a zone of transition from hght travel 
time dominated flashes to those dominated by diffusion. By 
setting tc — tsc and solving for one obtains an expression 
for the radius of a RSG progenitor which is in the transition 
zone (i?*Tz): 



R*TZ = 697 



Ei„ 



0.34cm2g 
0.68 



10^1 erg 



WMq 



Rg 



.0. 



(13) 



The duration of a flash from a transition zone supernova 
(tTz) is simply _R*tz/c: 



tTZ = 27.0 



0.34cm2g-i 



0.50 / \ 0.24 
Pi 



/ Sin 



V 1051 erg 



10M(7 







(14) 



slightly longer than the LOBSTER dwell time. 

Hence, for breakout bursts with duration t < txz, the 
progenitor radius is ct, and for those with longer duration, 
can be constrained via the tsc deflnition of MM 99. Sim- 
ilar equations can be derived to deflne the transition zone 
for BSG progenitors, though it should be noted that, for 
the parameters considered (Table the breakout flashes 
from BSGs are always light travel time dominated when 
i^in ^ 10^^ erg; hence, a transition zone exists only for BSG 
explosions with less than this canonical energy. 

As an example, consider a supernova with Ein = 10^^ 
erg which produces a breakout flash with a well-measured 
duration t. Using the parameters given in Table|5]and equa- 
tion 1141 one can establish minimum and maximum values 
for the light travel and diffusion times of RSG and BSG 
progenitors and draw conclusions from how t falls into these 
ranges. If t < 116 s, the progenitor is a BSG and its radius 
is ct (since no RSG progenitor has such a small light travel 
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Table 5. Constraints on progenitor type and _R* from timing 
alone. An explosion energy Ei^ = 10^^ erg is assumed. Note that 
1100 - 2080 s is the region spanned by the transition zone (TZ); 
hence, a breakout flash with this duration could be either tc or 
tsc dominated - more information is required. 



t (s) 


Progenitor 






type 




- 116 


BSC 


ct 


116 - 230 


need colour info 


ct 


230 - 1100 


RSG 


ct 


1100 - 2080 


RSG 


near TZ 


> 2080 


RSG 


can constrain using tse 



time). For a duration of 116 - 230 s, the colour of tiie flash 
could be used to distinguish between RSG and BSG (see 
il3.4L and the progenitor radius is ct (since both RSGs and 
BSGs can have a tc of this magnitude) . A flash > 230 s long 
denotes a RSG (since no BSG is large enough to produce 
a flash this long); for 230 s < t < 1100 s the radius is ct 
(since the minimum value of tTz is 1100 s), and for t > 2080 
s restrictions can be placed on i?* using the equation for t^c- 
Flashes with duration 1100 - 2080 s are RSG type and could 
be either tc or tsc dominated (since the maximum value of 
tTz is 2080 s) . Note, however, that the duration of a flash is 
not well constrained by LOBSTER if it exceeds the ~ 400 s 
dwell time of the experiment. 

These results are summarized in Tabled 



- ^ ^ X 1000 R 




N =0 1 2 3 5 7 9 

21 



-2 -1 (1 1 2 

10 10 10 10 10 

C3/C, 

Figure 4. Duration and C3—2 colour index for all progenitor 
models; the solid line indicates BSGs and the dashed line RSGs. 
Each curve corresponds to a value of TVjj, which increases to the 
right in units of lO^'^ cm~^ as indicated. Dotted lines mark con- 
stant radius. Models marked with "x" are uncharacterized and 
those with "o" are unobservable f >l3.2.1i . The discontinuity in 
slope in the RSG plot denotes the transition from light travel 
time to diffusion time in the burst duration. Note that lines of 
constant radius for t > ixz are horizontal only because M^j and 
Ein are held constant in each model at 15Mq and 10^^ erg, re- 
spectively. 



3.4 Constraints From Timing and Colour 

By considering the colour of a flash in addition to its dura- 
tion, more information about the supernova can be deduced 
and certain degeneracies may be broken. For instance, a flash 
with duration 116 s < t < 230 s (from a 10^^ erg supernova) 
is light travel time dominated, but both large BSGs and 
small RSGs produce flashes of such duration. The colour of 
the flash may help to break this degeneracy and allow for 
an identification of the progenitor type. 

Tse varies inversely with i?*; hence, larger progenitors 
produce redder flashes of lower radiation temperature with 
more channel 1 photons than smaller progenitors do. As can 
be seen in Figure |1] a RSG produces a flash which is dis- 
tinctly bluer in colour than that of a BSG of the same radius 
(same flash duration); however, RSG flashes in general are 
redder and brighter since RSGs typically have much larger 
radii than BSGs. The colour difference between BSG and 
RSG flashes of the same radius is enhanced by increasing 
the absorbing column density, as the redder BSG flashes 
lose proportionally more channel 3 counts than their RSG 
counterparts (this effect is evidenced by the divergence of 
the two types in C3-2 space). Figure|^shows that RSG and 
BSG flashes lose proportionally the same amount of counts 
in channels 1 and 2 as absorption increases (no divergence). 
It seems unlikely that the "type degeneracy" mentioned in 
>l3.3l can be broken without an estimate for A''h and an ac- 
curate measure of the flash colour. 

As previously noted, increased absorption causes the 
values of C2-1 and C3-2 to increase as lower energy pho- 
tons are preferentially absorbed. A greater Ein also leads 
to increased colour values due to its proportionality with 



1000 R 



1 \ ' . 




Figure 5. Duration and C2-1 colour index for all progenitor 
models; the solid line indicates BSGs and the dashed line RSGs. 
Each curve corresponds to a value of A^h, which increases to 
the right in units of lO^'^ cm~^ as indicated. Dotted lines mark 
constant radius. Points marked with "x" indicate models which 
are uncharacterized f il3.2.1i . The discontinuity in curvature of 
the RSG plot denotes the transition from tc-dominated to tgo- 
dominated duration. Note that lines of constant radius for t > txz 
are horizontal only because Mcj and E^^ are held constant in each 
model at 15Mq and 10^-'^ erg, respectively. 
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Figure 6. Colour-colour plot for all progenitor models; the solid 
line indicates BSGs and the dashed line RSGs. Each solid or 
dashed curve corresponds to a value of Nn , which increases to the 
right in units of lO^'^ cm~^ as indicated. Dash-dotted lines mark 
constant colour temperature Tbc for BSG models, dotted lines 
mark constant Tge for RSG models. With each Tse is provided (in 
parentheses) the value of i?^ at which this value is attained in the 
case Bin = lO^-' erg and M^j = ISA/q. Points marked with "x" 
indicate models which are uncharacterized ( i|3.2.1l . 



Tso, whereas Mej is inversely proportional to the radiation 
temperature (c.f. Table 0. The effect of Ein on colour is 
very small (~ 0.2 power-law dependence) compared to that 
of A'^H (exponential dependence), and the effect of Moj is 
minute (~ —0.05 power law dependence); -Ein and Mej do 
have a significant effect on tsc- In the case of bursts with 
t > tdwcii, whose durations cannot be determined by LOB- 
STER , the SN properties are best constrained using flash 
luminosity and colour ( il3.51 . 

Accurate observations of a flash's duration and colour 
can pinpoint the location of the flash on a t vs. C3-2 or t vs. 
C2-1 plot (Fig.2]or|3l and hence constrain the value of A'^h 
in addition to R. If A^h is known from other observations, it 
can be used in tandem with LOBSTER colour observations 
to constrain Mej and Ein- 

Potentially better constraints come from a colour-colour 
diagram. Figure |S] shows that the bluening effect of extinc- 
tion can, to some extent, be disentangled from the intrinsic 
flash colour to provide an estimate of both A^h and Tsc. The 
latter can then be taken as a constraint on _Ein, Mcj , and ii*, 
which is signiflcantly more illuminating if is constrained 
by other means. If, in addition, Esc/t is determined by com- 
paring the observed count rate to a value of D (from optical 
foUowup; Figs. |7| and |H| >I3.5II . then all explosion parameters 
are constrained. 



3.5 Constraints from Luminosity and Colour 

Flashes that outlast the LOBSTER dwell time have dura- 
tions that can only be constrained rather than measured; 
this puts a lower limit of roughly 170 Rq on the progenitor. 
Their durations may be either diffusion or light travel-time 




Figure 7. A theoretical X-ray Hertzsprung-Russell diagram of 
breakout flashes, in which luminosity (measured via the rate of 
observable photon production) and colour (here 03-2) are com- 
pared for RSG models. These quantities are observable even if 
the burst is not viewed for its entire duration, and thus are useful 
for flashes with t > tdwcll- Each curve corresponds to a value of 
A^H, which increases in units of lO'^^ cm~^ as indicated. Dotted 
lines mark constant values of Tge as indicated; Tge increases to 
the right along each curve. All models are evaluated at constant 
-Esc/t (10''^ erg/s); a change in -Esc/t results in a vertical shift 
of the curves (the dashed arrows shows a 0.9 dex shift). Mod- 
els marked with "x" are uncharacterized and those with "o" are 
unobservable (" i|3.2.1l . 



dominated (JU, but equation ll^Jl is unlikely to discriminate 
between these possibilities with only a lower bound on tse. 

The readily observable quantities for these long bursts 
are their X-ray colours and brightnesses. We will as- 
sume that optical follow-up allows a determination of 
flash distances: in that case, one can construct an X-ray 
Hertzsprung-Russell (H-R) diagram for them (Figure |7J. In 
the H-R diagram, the rate of observable photon production 
[Count rate]) is plotted versus a colour index. Spec- 
ifying Tse, -Esc/t, and Nn produces a point on the plot. We 
hold -Bse/t constant and vary Tsc and Nh to illustrate trends; 
varying Esc/t simply results in a vertical shift of the curves. 

By placing a point on Figure|7|representing the observed 
characteristics of a flash, it possible to flx Tsc- With other 
information, such as that provided by the colour-colour di- 
agram (Fig.|HJ, Nu and Ese/t can be constrained. Through 
the scaling equations Q, M^, Ein, and /?* can, in turn, 
be constrained. Although all BSG progenitors are expected 
to produce flashes with t < tdwcii, we include an H-R di- 
agram for BSGs for the purposes of comparison (Fig. (HJ. 
For higher Tsc models, BSG flashes are effectively indistin- 
guishable from RSG flashes of the same flash luminosity, 
-Esc/t. The H-R diagram does not present a good means of 
constraining progenitor type for flashes with t < tdwoii; the 
other methods described in this section work best. For long 
flashes which are not observed in their entirety, the H-R 
diagram could prove invaluable. 
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Figure 8. As in fig. |7| but for BSG models (tliick solid lines). 
RSG models appear as thin dashed lines for comparison. These 
quantities are observable even if the burst is not viewed for its 
entire duration, although all BSG progenitors are expected to 
have duration t < tdwell. Each curve corresponds to a value of 
A'^Hi which increases in units of 10^^ cm~^ as indicated. Dotted 
lines mark constant values of Tgo as indicated; Tse increases to the 
right along each curve. All models are evaluated at constant Esc/t 
(10*^ erg/s); a change in Ese/t results in a vertical shift of the 
curves (the dashed arrows shows a 0.9 dex shift). Models marked 
with "x" are uncharacterized and those with "o" are unobservable 



4 SHOCK TRAVEL TIME 

The travel time for the Shockwave responsible for the break- 
out flash, or equivalently, the time lag between the emer- 
gence of neutrinos and gravity waves from the exploding star 
and the emergence of shock breakout, is well approximated 
by 



dr 



= 1.26- 



E 



1/2 



/' 


m{r) 


1/2 


■po(r)r^" 




[ Mej J 




m(r) 



(15) 
dr 



Stellar mass and density profiles are needed to define Vs 
(Eq. EJ; the authors used progenitor structures provided by 
Woosley and Nomoto. For a IQMq, ASRq BSG with M<,j = 
14M0 and Sin = lO'^^ erg, At^,^ is 2 hours. For an ISA^q, 
4507?Q RSG with M^j = 16.5Mq and Ein = 1.1 x 10^^ erg, 
the delay between neutrino/gravity wave emission and shock 
breakout is 36 hours. It should be noted that this method to 
determine the time of shock emergence using equation I15II 
is considerably more accurate tha n the Sedov solution for a 
constant density envelope, used bv lWooslev et al.l i2002h and 
others. The determination of At^^-y is not of consequence ob- 
servationally unless the supernova in question occurs within 
a few kpc; in this case, it is feasible to detect the neu- 
trino/gravity wave emission with the appropriate detectors. 
Near-future gravity wave de tectors will only be able to de- 
tect SNe within our Galaxy dOtt et alJl2004) . whereas neu- 



trino signals can be detected out to the Magellanic clouds. 
Neutrino and/or gravity wave detection would provide: 1. 
early warning that a breakout fiash (and supernova) will oc- 
cur, and 2. a measurement of At^^.y, which could be used 
to further constrain the supernova properties. It should be 
noted, however, that the rate for Galactic supernovae is very 
low (~ 0.01 yr~^), that only ~ 9% of these fall within the 
LOBSTER field of view, and that most of these occur in 
the Galactic disk, where A''h > 10^^ cm~^ typically; hence, 
constraining SNe properties using shock travel time could 
only be done under very fortuitous circumstances. 



5 SHOCK BREAKOUT OBSERVATION WITH 
OTHER INSTRUMENTS 

In examining the capabilities of an instrument to observe 
the shock breakout flash, this paper has focused on the pro- 
posed LOBSTER wide-field X-ray detector; however, break- 
out hashes are also detectable by other current instruments 
- though less frequently and in less detail. Preliminary work 
by the authors indicates that the EPIC-PN instrument of 
the XMM-Newton X-ray space observatory, with a narrower 
field of view but much higher effective area than LOBSTER, 
is capable of detecting fiashes at moderate redshifts. The 
~ 0.2 deg^ field of view increases Dmin ( ^3. Ill , or equiva- 
lently, Zmin, to Zmin = 0.2; but the high effective area allows 
detections of the canonical fiashes to Zmax ~ 0.9 (c.f. Dmax, 
H3.H . The Chandra X-ray observatory's HRC-I has a simi- 
lar field of view, and hence Zmin, as XMM Newton, with a 
lower effective area and correspondingly lower «max ~ 0.2. 
The PSPC instrument on the 7?OS'j4r satellite, with a large 
field of view and intermediate effective area, has Zmin = 0.1, 
Zmax ~ 0.6. The assumed columns are A^h = 2 x lO'^^ cm~^ 
in the host galaxy and A'h = 0.5 x 10^^ cm~^ locally. These 
results indicate that there may be a number of undiscovered 
breakout flashes in the backgrounds of archived XMM New- 
ton and ROSAT images, of order 10 for each instrument. In 
the HRC-I images, there may be of order 1 BSG breakout 
flash and l i kely n one from RSGs; it should be noted that 
iRace et al.l (|2to3) have come to a similar conclusion with 
regard to the ACIS instrument on Chandra. 

Because the breakout flash durations are much shorter 
than the typical integration times for these instruments, and 
because a typical flash detection would be dimmer than the 
X-ray sources observed, the probability of a serendipitous 
discovery is almost zero. This may explain why the archived 
breakout flashes have gone unnoticed, even if present in the 
data. An in-depth search of the archived images would be 
required to locate the detected flashes. 

A sensitive optical survey with a sufficiently large field 
of view may be able to detect the Rayleigh- Jeans tail of 
the breakout flash spectrum. Preliminary calculations indi- 
cate that, to an extinctionless limiting magnitude of ~ 23 
mag in both the U and B bands (Johnson-Morgan system), 
there is one flash per square degree per 200/^^ years for the 
canonical RSG, and one flash per square degree per 600f~y^^ 
years for the canonical BSG. Extinction in the host galaxy 
typically causes the limiting magnitude to increase to ~ 24 
mag. 

We have ruled out the possibility of detecting high red- 
shift breakout flashes with optical telescopes; the high lu- 
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minosity distances involved more than compensate for the 
cosmological K-correction. 



6 DISCUSSION 

We have combined an analytical theory for the depen- 
dence of shock breakout parameters (MM 99) with numeri- 
cal simulations of shock breakout flashes felein & Chevalier 
| l978l: iBlinnikov et al.l Il998l: lEnsman fc Burrow j Il992t 
lBhnnito>^^nT200(lir to predict the expected signal ob- 
served by the LOBSTER spaceborne soft X-ray camera. Our 
emphasis has been on the reconstruction of supernova pa- 
rameters - primarily radius, mass, explosion energy, and ob- 
scuring column - from the data 

Supernova radius is the most tightly constrained of the 
three quantities f >13.3t . For all events which have a dura- 
tion that can be measured individually, because it is shorter 
than the dwell time of the instrument, the duration is set 
by the light travel time of the star rather than the leakage 
of photons from the surface layers of the explosion. These 
events come from stars with R ^ 170Rq, i.e., blue super- 
giants and relatively compact red supergiant progenitors. 
Such flashes are characteristically harder and dimmer than 
those from more extended RSG progenitors; nevertheless, 
they are visible to distances of several hundred Mpc and 
should be observed in the hundreds per year (Table |5J. 

For these stars, the flash colour temperature Tsc and 
absorbing column A^h can be estimated for extinctions 
A'h ^ 2 X 10^^ cm~^ by placing them on duration-colour 
(Figs. 131 and [SJ and colour-colour (Fig. IBJ diagrams. Tsc 
provides a constraint on the explosion energy and ejected 
mass through the equations derived by MM 99 (Table 0; 
however, these quantities cannot be derived independently 
without additional information. In the case that optical fol- 
lowup provides a distance, this degeneracy can be broken 
(Figure 0. 

Events that outlast LOBSTER'S dwell time are overrep- 
resented because they can be detected even if they do not 
begin while inside the fleld of view ( tl3.1ll . However they are 
also poorly characterized, because only a lower limit is avail- 
able for their duration. With optical follow-up they can still 
be placed on FigureEl allowing an estimate of the obscuring 
column and the colour and luminosity of the breakout flash. 



6.1 Implications 

From the above analysis we draw several conclusions. The 
LOBSTER space observatory will provide a census of Type 
II supernova events that is complete within ~ 250 Mpc (Ta- 
bleland FigureEl and up to columns Nh > lO^^cm^^. The 
number of events expected in three years of LOBSTER ob- 
servations is quite uncertain, but for a conservative estimate 
of host-galaxy extinction (Table 2J it ranges from ~ 50 to 
~ 600 depending on the typical radius of the progenitor 
stars. 

A population of dim Type II supernovae from blue 
progenitor s is suggested by SN 1987A and its analogues 
llSchmitz fc Gaskel 1988; F iflipenkdll983 lYoung fc BranchI 
[^8^ ~ ls diarfeij_^9^|)^_ although its sign ificance is contro- 
versial ll^^e^e^hTMcOuri [l9^. SN 1987A ana- 
logues should appear in the data as flashes whose duration is 



shorter than the dwell time of the instrument. These events, 
whose properties are difficult to infer from optical observa- 
tions due to the contribution of ^''Co ( tll.H . will provide 
strong constraints on stellar evolution at the point of core 
collapse. 

6.2 Caveats 

The current work relies on the extrapolation of breakout 
fiash properties for a variety of stellar progenitors from only 
a pair of numerical calculations, by means of the analyti- 
cal scaling relations given by MM 99; see il3.2l It is possible 
that elements of the radiation dynamics cause the overall 
shape of the breakout spectrum to change across this range, 
which would introduce a systematic error into our predic- 
tions of X-ray observables. This can only be tested with a 
more systematic survey of breakout properties, preferably 
using multigroup radiation hydrodynamics simulations. 

One limitation of the MM 99 scaling relations is that 
they cover only two possible forms of the outer density pro- 
file (p oc depth"): n = 3/2, representing convective en- 
velopes, and n — 3, representing radiative envelopes with 
constant opacity. Other possibilities should be considered. 
For instance, convection is inefficient near the surfaces of 
red supergiants, and this effects shock propagation in the 
outermost layers (as MM 99 noted). A superadiabatic layer 
has a value of n that is lower than the adiabatic value. On 
the other hand, the MM 99 analysis and our own investiga- 
tions show that the detailed structure of the outer envelope 
plays a rather minor role in the breakout fiash intensity, 
colour, and duration (see also the Appendix). 

6.3 Speculation: Asymmetric Explosions 

As noted above, the short breakout flashes that fit within the 
LOBSTER dwell time have durations set by the star's light 
travel time. This makes it possible for asymmetries in those 
explosions to affect the time dependence of their breakout 
flashes. A spherically symmetric explosion will exhibit a pro- 
gression in brightness and colour that represents the growth 
of the emitting area and the change of its limb darkening 
as the observable portion of th e breakout moves from th e 
front to the side of the star fe.g.. lEnsman &c Burrow3ll992l) . 
If the explosion is sufficiently asymmetric, then this pattern 
will be disturbed by the motion of the shock front across the 
face of the star. Unfortunately, the most compelling evidence 
would derive from an observation of time-dependent linear 
polarization of the emerging X-rays; a difficult quantity to 
observe. 

Asymmetries can derive from several sources: asymme- 
tries in the envelope distribution (due to rotation and con- 
vective eddies); the growth of (weakly) unstable perturba- 
tions in accelerating shocks; and the asymmetry of the cen- 
tral engine driving the explosion, e.g., in the case of a jet- 
driven explosion. (Note that gam ma-ray burst je ts are not 
likely to escape supergiant stars iMatzneJliooai : however 
they may imprint an asymmetry on the explosions.) 

Concentrating on the first of these, how much asymme- 
try would be required to perturb the shock travel time by 
an amount comparable to the light travel time? The shock 
speed is roughly 1(10MQ/Mej)% of c on average, so a rela- 
tive difference in travel time of the same order is required; 
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this would arise from a comparable asymmetry in the pro- 
genitor. This degree of asymmetry may exist in blue super- 
giants if, for instance, they have undergone tidal interactions 
with companion stars in previous red supergiant phases or 
are close to their Roche radii at the time of explosion. 



ACKNOWLEDGEMENTS 

We thank the referee, Roger Chevalier, for pointing out that 
early SN luminosity provides a constraint on stellar radius. 
We are grateful to Nigel Bannister and the LOBSTER sci- 
ence team for their correspondence and for specifications of 
the instrument, and to Neil Brandt, David Ballantyne, and 
John Monnier for helpful suggestions. AJC was supported 
in part by an NSERC undergraduate fellowship. CDM is 
supported by NSERC and the Canadian Research Chairs 
Program. 



REFERENCES 

Bannister, N. 2003, private communication. 

Blinnikov, S. et al. 1998, ApJ, 496, 454. 

Blirmikov, S. et al. 2000, ApJ, 532, 1132. 

Blaauw, A. 1991, NATO ASIC Proc. 342: The Physics of 

Star Formation and Early Stellar Evolution, 125 
Cappellaro, E. et al. 1997, A&A, 322, 431. 
Cappellaro, E. et al. 1999, A&A, 351, 459. 
Chevalier, R. A. 1992, ApJ, 394, 599 

Colgate, S. A. 1968, Canadian Journal of Physics, 46, 476 

Colgate, S. A. & White, R. H. 1966, ApJ, 143, 626 

Ensman, L. & Burrows, A. 1992, ApJ, 393, 742. 

Fillipenko, A. V. 1988, Supernova 1987A in the Large Mag- 
ellanic Cloud, 106 

Fukugita, M. et al. 1998, ApJ, 503, 518. 

Gandel'man, G. M. and Frank-Kamenetsky, D. A. 1956, 
Sov. Phys. Dokl., 1, 223. 

Grassberg, E. et al. 1971, Ap&SS, 10, 28. 

Grover, R. & Hardy, J. W. 1966, ApJ, 143, 48 

Hamuy, M. 2003, ApJ, 582, 905. 

Iglesias, C. A. & Rogers, F. J. 1996, ApJ, 464, 943 

Imshennik, V. S. and Nadezhin, D. K. 1989, Sov. Sci. Rev. 
E. Astrophys Space Phys, 8, 1. 

Klein, R. & Chevalier, R. 1978, ApJ, 223, L109. 

Klein, R. et al. 1979, ApJ, 234, 566. 

Litvinova, I. & Nadezhin, D. 1985, Soviet Astron. Lett., 11, 
145. 

Lundqvist, P. & Fransson, C. 1988, A&A, 192, 221 
Matzner, C. & McKee, C. 1999, ApJ, 510, 379. 
Matzner, C. D. 2003, MNRAS, 345, 575 
Ott, C. et al. 2004, ApJ, 600, 834. 
Popov, D. 1993, ApJ, 414, 712. 
Priedhorsky, W. et al. 1996, MNRAS, 279, 733. 
Race, D. et al. 2003, AAS Meeting 203, 142.02. 
Sakurai, A. 1960, Comm. Pure Appl. Math., 13, 353. 
Schaefer, B. E. 1996, ApJ, 464, 404 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 
500, 525 

Schmitz, M. F. & Gaskel, C. M. 1988, Supernova 1987A in 
the Large Magellanic Cloud, 112 
Spergel, D. N. et.al. 2003, ApJS, 148, 175. 



Stritzinger, M. et al. 2002, AJ, 124, 2100 
Tan, J. C, Matzner, C. D., & McKee, C. F. 2001, ApJ, 
551, 946 

van den Bergh, S. 1991, Proceedings of the Astronomical 

Society of Australia, 9, 13 
van den Bergh, S. & McClure, R. D. 1989, ApJ, 347, L29 
Van Dyk, S. et al. 2003, PASP, 803, 1. 
Wilms et al. 2000, ApJ, 542, 914. 
Woosley et al. 2002, Rev. Mod. Phys., 74, 1015. 
Young, T. R. & Branch, D. 1988, Nat, 333, 305 



APPENDIX A: STELLAR OUTER DENSITY 
COEFFICIENTS 

The parameter pi /p* , which describes the density structure 
of the outermost regions of a progenitor star O, appears 
in the MM 99 equations for the properties of the breakout 
flash. Though the near-surface structure plays little role in 
determining the temperature and energy of a breakout flash, 
it is significant in dictating the duration tsc of a diffusion- 
dominated fiash. 

The following sections we estimate pi/p* for the ra- 
diative envelopes of blue supergiants and for the convective 
envelopes of red supergiants. In each case we motivate re- 
placing pi / p* with a particular numerical value; only in ex- 
ceptionally well-observed breakout fiashes could the depen- 
dence of pi/p* on stellar parameters be used to reconstruct 
them. 



Al Outer Density Coefficients for Blue 
Supergiants 

It is possible to re-write the MM 99 scaling equations for 
BSGs (Table such that pi/p* is eliminated in favour of 
the mass, luminosity, and compositional parameters of the 
progenitor star. For progenitors with a radiative envelope 
(i.e., BSGs), pi may be expressed as 



pi _ aifiniH)'^ G'Hd't 
P* 



192fc| 



A/ej 1 



(Al) 



where 1 — /3 = I/*/I/Edd is the ratio of the stellar luminosity 
to the Eddington limit. For Thompson opacity. 



1 - /3 = 0.259 



which leads to 
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IOMq ^ 
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Using theoretical stellar model data from various au- 
thors, we have established a semi-empirical mass-luminosity 
(M-L) relation for presupernova stars. This M-L relation, 
show in Figure IXTl is 
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Table Al. Alternate scaling equations for progenitors with 
radiative envelopes, without using the M-L relation. These 
should be read as exponents in formulae, e.g., tae = 



101-^5[/k/(0.34 cm^g-i)] 



-0.14 



cgs units. 



Parameter 




Esc 




10 


6.11 


47.37 


1.65 


/k/0.34 cm^ g-l 


-0.14 


-0.84 


-0.45 


An/lO^^ erg 


0.18 


0.58 


-0.72 


Mej/1OM0 


-0.11 


-0.36 


0.46 


i?*/5OR0 


-0.48 


1.68 


1.90 


M*/1OM0 


0.18 


-0.22 


-0.73 


1 — L*/iEdd 


0.18 


-0.22 


-0.73 


/i/0.62 


0.18 


-0.22 


-0.73 


il + XH)/1.7 


-0.046 


0.054 


0.18 
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Figure Al. Luminosity- mass relation for presupernova stars. 
Plotted are the terminal masses and luminosities for the Padova 
tracks with solar-scaled metallicity Z = 0.008 (Salachnich et al. 
2000, A & A 361, 1023, crosses), for Geneva tracks of the same 
metallicity (Scherer et al. 1992 A & AS Supp. 98, 523), and 
for a progenitor model for SN 1987A (Hashimoto, Shigeyama, 
& Nomoto 1989, circle). Also plotted is a simple power law fit 
fEa.lX4t. 
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and it is valid for the mass range 10 < M-^/Mq < 20. Sub- 
stituting the M-L relation into equation IIA3II results in: 



Table A2. Alternate scaling equations for progenitors with ra- 
diative envelopes, using the mass-luminosity relation (Eq. IA4I '). 
These should be read as exponents in formulae, e,g., tse = 



lO^'''* [/K/(0.34cm2g-i)]- 



cgs units. 



Parameter 


^SC 


Esc 




10 


6.08 


47.41 


1.78 


/k/0.34 cm2 g-l 


-0.14 


-0.84 


-0.45 


Bin/lO^l erg 


0.18 


0.58 


-0.72 


Mej/1OM0 


-0.11 


-0.36 


0.46 




-0.48 


1.68 


1.90 


M*/1OM0 


0.10 


-0.12 


-0.40 


1 — L*/ L^dd 


0.18 


-0.22 


-0.73 


m/0.62 


0.18 
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-0.73 


(1 + Xh)/1.7 


-0.046 
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0.18 
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From equation ljA5|l it follows that pi/p* varies almost 
linearly from ~ 0.13 to ~ 0.26 for 10 < PL/Mq < 20, 
assuming a 1.5Mq remnant and << Lsdd- Thus, we will 
adopt the average, ~ 0.2, as a fiducial value for all BSG 
progenitors. 



A2 Outer Density Coefficient for Red Supergiants 

Combining equations (9), (14), and (48) from MM 99 results 
in: 

> 5/2 



12.17(7) (0.37 - O.lSq + 0.096g 



(A6) 



^ = (12.62 

P* 

where q = 1 — Mcm/K'h and Mcnv is the mass of the outer 
stellar envelope. In this case, p\/p* exhibits an approxi- 
mately 1/q dependence for the most plausible range of q 
values, 0.3 < g < 0.5 (see Figure 5 of MM 99). For these 
values of q, pi/p* varies from ~ 0.54 to ~ 0.33; we adopt 
0.5 as the standard for RSGs. 



